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Abstract 

One of the most remarkable results of quantum mechanics is the fact that many-body quantum systems 
may exhibit phase transitions even at zero temperature [1J. Quantum fluctuations, deeply rooted in Heisen- 
berg's uncertainty principle, and not thermal fluctuations, drive the system from one phase to another. Typ- 
ically, the relative strength of two competing terms in the system's Hamiltonian is changed across a finite 
critical value. A well-known example is the Mott-Hubbard quantum phase transition from a superfluid to 
an insulating phase El |3l, which has been observed for weakly interacting bosonic atomic gases. However, 
for strongly interacting quantum systems confined to lower-dimensional geometry a novel type of quantum 
phase transition may be induced for which an arbitrarily weak perturbation to the Hamiltonian is sufficient 
to drive the transition HI [3. Here, for a one-dimensional (ID) quantum gas of bosonic caesium atoms 
with tunable interactions, we observe the commensurate-incommensurate quantum phase transition from 
a superfluid Luttinger liquid to a Mott-insulator|[6l|7]]. For sufficiently strong interactions, the transition is 
induced by adding an arbitrarily weak optical lattice commensurate with the atomic granularity, which leads 
to immediate pinning of the atoms. We map out the phase diagram and find that our measurements in the 
strongly interacting regime agree well with a quantum field description based on the exactly solvable sine- 
Gordon model[8]. We trace the phase boundary all the way to the weakly interacting regime where we find 
good agreement with the predictions of the ID Bose-Hubbard model. Our results open up the experimental 
study of quantum phase transitions, criticality, and transport phenomena beyond Hubbard-type models in 
the context of ultracold gases. 
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Ultracold atomic gases are a versatile tunable laboratory system for the investigation of com- 
plex many-body quantum phenomena[9|. The study of quantum phases and quantum phase tran- 
sitions is greatly enriched by the possibility to independently control the kinetic energy and the 
interactions. In deep optical lattice potentials the many-body dynamics for a weakly interacting 
gas is, to a very good approximation, governed microscopically by a Hubbard Hamiltonian[|2] 
with a local onsite interaction energy U and kinetic energy J, which corresponds to tunneling of 
atoms from one lattice site to the next. Experiments with Bose-Einstein condensates (BEC) of Rb 
atoms have demonstrated the quantum phase transition from a superfluid phase for large J to an 
insulating Mott-Hubbard (MH) phase[3]. The transition between these two phases was obtained 
by quenching 7 in a lattice of finite depth. Recent experiments with fermionic atoms have demon- 
strated the presence of a fermionic MH insulating state lfTOlfTTI . potentially opening the way to the 
study of high-temperature superconductivity in proximity of the MH phase in 2D. 

While the focus in the study of quantum phase transitions in the context of ultracold atoms 
has so far been on Hubbard-type physics in the weakly interacting regime, novel quantum phe- 
nomena occur in lower dimensions, where the effects of quantum fluctuations and correlations are 
enhanced. In a ID bosonic gas, strong repulsive interactions lead to the formation of a Tonks- 
Girardeau (TG) gas, where bosons minimize their interaction energy by avoiding spatial overlap 
and acquire fermionic properties [fT2] - iT5ll . The addition of an arbitrarily weak lattice potential 
commensurate with the atomic density, i.e. n ~ 2/ A, where n is the linear ID density and A/2 is 
the lattice periodicity, is expected to lead to a novel kind of quantum phase transitionJH [6j : the 
strongly correlated ID gas is immediately pinned by the lattice and the superfluid TG phase is 
turned into an insulating, gapped phase. Figure [T] contrasts the Hubbard-type superfluid-to-Mott- 
insulator transition to this pinning transition. Given the universality of ID quantum physics, the 
pinning transition will occur for interacting bosons as well as for fermions in ID and has been 
discussed with respect to a variety of quantum models in low dimensions [4J. 

The pinning transition is described by the (1+1) quantum sine-Gordon (sG) model, which is an 
exactly solvable quantum field theory, extensively studied in high energy, condensed matter, and 
mathematical physics[5J. The sG Hamiltonian reads 

^ = S / dx ^ dx ^ 2 + (<9 ^ )2 + ^ cos (v^e)]. (i) 

Here, d x 6 and d x (j) are the fluctuations of the long-wavelength density and phase fields 6 
and 0, respectively, of the hydrodynamic description of the ID liquid with commutation rela- 
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Bose-Hubbard Sine-Gordon 



FIG. 1 : Comparing two types of superfluid-to-Mott-insulator phase transitions in ID. Schematic den- 
sity distributions (grey) in the presence of a periodic potential (red solid line), a, Mott-Hubbard type quan- 
tum phase transition for weak interactions @. The system is still superfluid at finite lattice depth (top). The 
transition to the insulating state is induced by raising the lattice depth above a finite critical value (bottom), 
b, Sine-Gordon type quantum phase transition for strong interactions [6]. In the absence of any perturba- 
tion, the system is a strongly correlated superfluid (top). For sufficiently strong interactions, not necessarily 
infinitely strong, an arbitrarily weak perturbation by a lattice potential commensurate with the system's 
granularity induces the transition to the insulating Mott state (bottom). 

tion [d x 0(x), (j>(y)] = ilt5(x — y), v s is the velocity of the soundlike excitations of the ID gas, 
V = Vnnj (hv s ) is proportional to the depth V of a weak lattice^ [6j, and ft is Planck's constant ft 
divided by lit. For vanishing lattice Y = 0, Eq. ([T]) describes a Luttinger liquid, where the strength 
of interactions is parameterized by the dimensionless parameter K = ftnn / (mv s ) , which determines 
the long-distance power-law decay of the correlation functions, e.g. (n(x)n(x r )} ~ l/\x — x'\ 2K , 
with m the atomic mass. The sG model with a weak but finite lattice predicts a quantum phase tran- 
sition of the Berezinskii-Kosterlitz-Thouless (BKT) type between a superfluid state for K > K c = 2, 
where the shallow lattice is an irrelevant perturbation, to an insulating Mott phase for K < K c , for 
which the spectrum is gapped for any value of Y. 

While in general K is a phenomenological parameter, in the case of a ID bosonic gas it can be 
microscopically related to the Lieb-Liniger parameter y = mgj (ft 2 n), which characterizes interac- 
tions in a homogenous ID system lfT6l (see Methods). Here, g ~ 2ftco±a^ is the coupling constant 
of the ID 5-function interaction potential U (x) = g8(x), where co± is the frequency of transverse 
confinement and is the 3D scattering length. The strength of interactions, and thus K, can be 
tuned by varying ci^d near a Feshbach resonance lfTTll . The TG regime corresponds to y 3> 1. Using 
the relation between K and y, Biichler and coworkers [6 J have shown that particles are pinned for 
experimentally accessible values of y > y c ~ 3.5 in the limit of a vanishingly weak lattice. The 



4 



pinning transition is expected to continuously transform into the MH-type quantum phase transi- 
tion, which occurs for the weakly interacting gas when the lattice depth becomes sufficiently large. 
Here, using a quantum gas of caesium (Cs) atoms with tunable interactions confined to an array 
of independent ID tubes (see Methods), we drive the superfluid-to-Mott-insulator phase transition 
by varying y and determine the phase boundary all the way from the strongly to the weakly inter- 
acting regime using modulation spectroscopy and measurement of transport. For shallow lattices 
under conditions of commensurability, we observe immediate pinning of the particles for strong 
interactions when y > y c . 

We first discuss our experiments in the strongly interacting regime. We start with a 3D Bose- 
Einstein condensate (BEC) of typically 1.3 x 10 5 Cs atoms without detectable thermal fraction in 
a crossed-beam dipole trap with magnetic levitation[18J and initialize our system by creating a 
conventional 3D MH-state in a deep 3D lattice at U / (67) ~ 75 with precisely one atom per lattice 
siteflH. The array of ID tubes is obtained by reducing the lattice depth V along one direction. 
Our procedure ensures that a majority of tubes has a near-commensurate number density (see 
Methods). A Feshbach resonance allows us to control with a precision of 3 ciq limited by the 
presence of the magnetic field gradient. Here, ao is Bohr's radius. For the case of the shallow 
lattice, we probe the state of the system by amplitude modulation spectroscopy [[19]]. We determine 
the presence of an excitation gap E g by testing whether energy can be deposited into the ID system 
at a given excitation frequency /. The lattice depth V is modulated at / by 25% to 45% for 40 — 60 
ms. After ramping down the lattice beams adiabatically with respect to the lattice band structure 
and after a levitated expansion time of 40 — 60 ms[18], we detect the atoms by time- of- flight 
absorption imaging. We determine the spatial width of the atomic sample from a gaussian fit to 
the absorption profile and obtain the change 8 of the spatial width compared to the unmodulated 
case as a function of /. Two typical measurements are shown in Fig. |2fa), one in the superfluid 
phase and one deep in the ID Mott phase at the same value for the lattice depth, V = 1.5(1)75r, 
where Er = h 2 / (2mA 2 ) is the photon recoil energy. For weak interactions the system exhibits a 
linear increase for 8 as a function of /, which we attribute to the superfluid character of the gas. 
For strong interactions, the increase, after a slow rise, shows a clear kink. We attribute the initial 
slow rise to excitation of residual superfluid portions of our inhomogeneous system and the sudden 
change in slope to the presence of an excitation gap. We associate the axis intercept f g obtained 
from a linear fit to the steep part of the spectrum with the frequency of the gap. To determine 
the phase transition from the ID Mott state to the superfluid state, we repeat this measurement for 
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FIG. 2: Modulation spectroscopy on bosons in ID. a, b, d, Excitation spectra for low, intermediate, 
and high lattice depth V. The change 8 of the spatial width after amplitude modulation is plotted as a 
function of the modulation frequency / for different values of 7 a, Characteristic spectra for V = 1.5(1)£r, 
one in the superfluid regime (squares, = 115(2) ao, y = 1.0(1)) and one in the Mott regime (circles, 
= 261(2) ao, 7 = 3.1(2)). The solid lines are linear fits to the high-frequency part of the spectrum. 
We determine the axis intercept f g as indicated, b, Spectra for V = 3.0(2)£r. The system is superfluid at 
7 = 0.51(6) (squares), while it exhibits a gap for y = 1.6(1) (triangles) and 7 = 4.1(3) (circles). Additional 
structure appears at higher frequencies due to the presence of band structure, c, Determination of the 
transition point for the case of the shallow lattice with V = 1.5(1)£r. The frequency f g is plotted as a 
function of y. The transition point is identified with the abrupt step. The solid line is an error-function 
fit to the data from which we determine the critical value y c y. The inset plots f g as a function of y for 
V = 3.0(2)£r. We identify the transition with the onset of the linear rise in f g . d, Spectra for V = 9.0(5)£r. 
Here, we give / in units of U. For comparatively weak interactions (squares, y = 0.10(3)) the spectrum 
exhibits abroad distribution characteristic of a superfluid (SF). For stronger interactions (circles, y = 8.1(4)) 
pronounced peaks develop, indicating the presence of a Mott insulator (MI). In a, b, d, the error bars for 8 
reflect the la statistical error. In c, the error bars for f g are derived from the la error on the fit parameters. 
The error for 7 results from the la statistical error of the independent input variables and the spread of 7 
due to the distribution of tubes. 6 




FIG. 3: Transport measurements on the ID Bose gas. Center-of-mass displacement xq as a function of am 
for different values of V (V = 9.0(5)£ R (diamonds), V = 5.0(3)£ R (squares), V = 2.0(1)£ R (circles)). We 
extrapolate the linear slope at small values for a^D and associate the transition point with the axis intercept. 
For the data with V = 2.0(1 )£r transport is not fully quenched as the condition of commensurability is 
not fulfilled for all atoms. All errors are the la statistical error. The inset plots the mesured critical ratio 
(U /J) c at the transition point as a function of lattice depth V. The dashed line indicates the theoretical result 
(U/J) c « 3.85 for the ID Bose-Hubbard regimeflB]]. 

a given depth V as we scan 7 by changing ci^d- A typical result is shown in Fig. [2jc). The gap 
closes as 7 is reduced. For values V < 2.0£r, we determine the critical value y c y at which the 
transition happens by an error-function fit to the data. Note that we always observe some small 
residual value for f a of about 120 Hz for weak interactions. 

For comparison, we present in Fig. ^b) and (d) excitation spectra for an intermediate value of 
the lattice depth and for the case of a deep lattice, respectively. For V = 3.0(2)£r the spectrum 
shows additional structure for high frequencies as band structure comes into play. We find that for 
V > 2.0£r the gap opens up approximately linearly as a function of 7 beyond a critical y c y, see 
inset to Fig.[2jc). For deep lattices we recover the discrete excitation spectrum of the Mott phase 
in the Hubbard regime ll3l [T9ll with a pronounced peak at / = 1.0 U/h. Additional peaks[20J can 
be found at / = 0.5 U/h and above f= 1.5 U/h. 
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For the case of a deep lattice, we find that the state of the system is very sensitively probed 
by transport measurements I12T1 l22l . A characteristic property of the Mott state is the inhibition 
of particle motion. In our experiment with the capability to tune interactions we expect the phase 
transition to manifests itself, at fixed V, through a strong suppression of transport when the strength 
of the interaction is raised above a certain critical value. Essentially, we test whether momentum 
can be imparted to the ID system as a function of interaction strength. For a given V we apply a 
weak axial magnetic force for a brief time to the interacting system, chosen such that the imparted 
momentum would be approximately 0.2M if the system were non-interacting. Then, as a function 
of a3D, we determine the center-of-mass displacement xq of the sample after a fixed time of flight. 
Fig. [3] shows that xq decreases monotonically with a^. For the case of a deep lattice with V = 
9.0(5)£r the quenching of transport is abrupt. At a certain critical value for a^D transport is fully 
inhibited li23l l24l . We find the critical by a linear fit to the decreasing data and by determining 
the axis intercept and derive from this a critical y c y ■ Reducing the lattice depth to V = 5.0(3)£r 
and V = 2.0(1 )£r leads to a less abrupt quenching of transport. For stronger interactions, the 
decrease starts to level off. Nevertheless, the initial decrease is still linear, allowing us to determine 
the critical y c .y by an extrapolation of the initially linear decrease to zero. The inset to Fig.[3]shows 
the measured critical ratio (U /J) c determined by our transport method as a function of lattice depth 
V. When we compare our results with the predicted value 11251 of (U /J) c ~ 3.85 for the transition 
in ID, we find a slight systematic overestimation of the transition point. This, however, is expected 
in view of e.g. the spatial inhomogeneity of the sample and the BKT-type nature of the transition 
in a finite size system. 

We summarize our results in Fig. [4] where we present the phase diagram as a function of 1 jy 
and V. The set {y c ,y} defines the phase boundary between the ID Mott insulator and the ID super- 
fluid. The measurements based on modulation spectroscopy cover a range from V = 4Er down to 
0.5£r (circles), while the transport measurements extend from V = 2£r to 10£r (squares). In the 
weakly interacting regime, l/y > 2, our data are in good agreement with the prediction of the MH 
model (dashed line). In the strongly interacting regime, l/y < 1, the measured phase boundary 
extrapolates to a finite critical value 1 /y c for the Lieb-Liniger parameter as the lattice depth V is 
reduced to zero. Our results are in excellent quantitative agreement with the theory for a commen- 
surate system based on the sine-Gordon model (solid line, see Methods), for which y c = 3.5. We 
also find good agreement between our two types of measurement techniques in the intermediate 
regime (V = 2£r to 4£r). Our results demonstrate the striking consequence of strong interactions 
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FIG. 4: Phase diagram for the strongly interacting ID Bose gas. Superfluid and Mott insulating phases in 
ID versus inverse Lieb-Lininger interaction parameter 1/y and optical lattice depth V in units of the photon 
recoil energy £r. The critical interaction parameter is y c . The data points (squares and circles) indicate 
the measured phase boundary. For strong interactions and shallow lattices we determine the transition by 
amplitude modulation spectroscopy (circles). Our measurements are in good agreement with the theory 
based on the sine-Gordon model (solid line, see Methods). For y > y c an arbitrarily weak lattice induces 
the transition to the insulating state. For weak interactions and deep lattices we probe the phase boundary 
by transport measurements (squares). Here, our measurements agree well with the prediction from the ID 
Bose-Hubbard model (dashed line). For the modulation measurements, the error in y is derived from the 
la error of the fit parameters. For the transport measurements, the error in y results from the la statistical 
error of the independent input variables and the spread of y due to the distribution of tubes. The inset plots 
the measured gap energy E g = hf g as a function of V for y = 1 1 ( 1) and compares our data to the analytical 
result for finite y as given by the sine-Gordon model (solid line, see Methods). Also shown is the universal 
behavior E g = V/2, which is valid for non-interacting fermions (dashed line). 
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in ID geometry in the presence of a lattice: Beyond a critical value y c , an insulating Mott state 
exists for vanishingly small lattice depth V. The particles are immediately pinned by the lattice. 

We measure a finite gap energy E g for 7 > y c in the regime of a shallow lattice. In the limit 
of 7 — > 00 and V — > one would expect the simple relation E g = V /2 as the bosonic system has 
become fully fermionized and the lattice effectively induces a band insulator of fermions[|6]|. In the 
inset to Fig. [4] we plot the measured E g as a function of V at fixed 7 = 11(1). For V < 1£r our data 
is in good agreement with the analytical result for the gap energy at finite 7 (see Methods). Note 
that, for V > 1£r, we observe a deviation for E g away from the predicted values. This deviation 
occurs at rather shallow lattices. However, one does expect the curve to have a reduced slope for 
deeper lattices, for which E g becomes of order U and is only weakly dependent on V. 

Our results are a benchmark realization of quantum field theory models with tunable parameters 
in cold atomic systems. These results open up the experimental study of the out-of-equilibrium 
properties of sine-Gordon-type models. In particular, thermalization in integrable models beyond 
the Luttinger liquid model, quenches across quantum phase transitions, and their relations to the 
breakdown of the adiabatic theorem in low dimensions can now be investigated with full tunability 
of system parameters. 
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METHODS 



ID Bose gas in a weak optical lattice. 



In the absence of the optical lattice, V = 0, the Luttinger liquid parameter K can be expressed 
in terms of the Lieb-Liniger parameter 7 = gm/ifp'n) for all strengths of interactions lfT6l l26l . 



For 7 < 10 and 7 > 10 one gets K ~ n/\Jy— y i l 2 /{2n) and K ~ (1 + 2/y) z , respectively. The 
addition of a weak but finite commensurate optical lattice with V < IE R realizes the effective 
sine-Gordon Hamiltonian Eq. ([T]). Using a perturbative renormalization group approach, the BKT 
transition line between the superfluid and the Mott-insulating phases can be derived in terms of V 
and 7 = 7c,v as 

V J - I 

£r Vv / r-r 3/2 /(2i) 

For small lattice depths, the integrable structure of the sine-Gordon model 11271 12811 allows one to 
derive the following analytical expression for the dependence of the spectral gap E g on V and K 



«T\ 



kK 



2(2-/0 



(k 2 v\ r[i-f] 



1 

2-X 



Vi6£ R ; r[i + f] 

1 , the dependence of the gap on V is 



Here, T is the gamma function. For strong interactions K 
linear, and E g approaches the free fermion value E g = V/2. In the vicinity of K = 2, the gap closes 
exponentially approaching the BKT transition line. 



Deep lattice: the Bose-Hubbard model. 

In the weakly interacting regime 7 <C 1, for V ^> 1£r, when all atoms occupy the lowest vi- 
brational state in each potential well of the lattice, the system can be described by the following 
Bose-Hubbard model[2J 

H = -J^(bp i+1 + h.c.) + ^b]bjb l b l . 

i i 

Here, Z?,- (bj) is the operator destroying (creating) a bosonic particle at the position of the z' th - 
well, J = 4£ R ( V /E R ) i exp [- 2 <J V /E R ] / 0r is the hopping energy, and U = V2ng(V /£r) 1/4 / X 
is onsite interaction energy. The quantum phase transition between a superfluid and a MH state 
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occurs at ll2"5l (U/J) c ~ 3.85, which determines a transition line in the (V, 7) - plane via 
Lattice loading and array of ID tubes. 

We create a 3D optical lattice by interference of 3 pairs of counterpropagating dipole trap laser 
beams at wavelength X = 1064.5 nm with l/e 2 beam waists of ~ 350 /im. The atomic BEC, 
initially trapped in a crossed-beam dipole trap, is adiabatically transferred to the 3D lattice by 
exponentially ramping up the power in the lattice laser beams within 300 ms. We create a 3D 
Hubbard-type Mott insulator with precisely one atom per site in the central region of the trap by 
adjusting the external dipole trap confinement. The array of vertically oriented tubes is created by 
ramping down the power in the vertically propagating beam pair with about 60 atoms in the center 
tube. Typical trapping frequencies for the tubes are 0) rz = 2n x (12300(200), 21.9(3)) Hz along 
the transversal and longitudinal directions, respectively. In view of our inhomogeneous system we 
calculate 7, for a given tube, by assuming a ID Thomas-Fermi distribution and taking the center 
density. The reported 7 is a weighted average over all tubes. We calibrate the depth of the lattice 
along the tubes by the pulsed Raman-Nath technique [|29l . 

Commensurability. 

To observe the pinning transition it is not necessary to fulfill the condition of commensurability 
precisely O. A finite commensurability parameter Q = 2%{n — n c ) corresponds to a shift 5/1 of the 
chemical potential. Here, n c = 2/X is the commensurate ID density. The system stays locked to 
the MI phase as long as 5/1 remains smaller than the energy necessary to add another atom. When 
Q rises beyond a critical value <2c(7, V), the system develops finite density excitations, which de- 
stroy the long range order of the MI. We find that, for the array of ID tubes, the commensurability 
condition in the superfluid regime is fulfilled best when the total atom number is chosen in such a 
way that the peak density of the center tube is approximately 1.2 n c . 
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